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We characterize the geometrical and topological aspects 
of a dynamical system by associating a geometric phase with 
a phase space trajectory. Using the example of a nonlinear 
driven damped oscillator, we show that this phase is resilient 
to fluctuations, responds to all bifurcations in the system, 
and also finds new geometric transitions. Enriching the phase 
space description is a novel phenomenon of "geometrical local- 
ization" which manifests itself as a significant deviation from 
planar dynamics over a short time interval. 
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Anholonomy is a phenomenon that arises if a quantity 
fails to recover its original value, when the parameters on 
which it depends are varied round a closed path. First 
introduced by Berry [1] in quantum mechanics as a geo- 
metric phase associated with the wave function, it is now 
recognized that this phase can emerge in classical sys- 
tems as well. The universality underlying this concept 
has proven to be very useful in a variety of systems [2,3]. 

In this Letter we present, for the first time, a formal- 
ism which uses this concept to quantify the geometric 
characteristics of complex phase space trajectories of sys- 
tems exhibiting chaotic dynamics. Here, the root of the 
anholonomy lies in the non-planarity of the trajectories. 
Using the example of a nonlinear driven damped oscilla- 
tor, we demonstrate that the geometric phase is a robust 
order parameter for the system, which not only responds 
to all bifurcations but also finds new geometrical tran- 
sitions. These include geometrical resonances leading to 
both discontinuous as well as continuous changes in the 
geometry of periodic and strange attractors. Our formal- 
ism integrates phase space dynamics and its geometrical 
aspects by including information about the non-planarity 
of the trajectory. Underlying this enrichment is the novel 
phenomenon of geometrical localization: the temporal lo- 
calization of the geometrical variables, when the devia- 
tion from planar dynamics occurs over small time inter- 
vals. 

Although our formalism can be generalized to higher 
dimensions, we consider (as an illustrative example) a 
dynamical system described by three first-order coupled 
differential equations, 



r t = v(r, a) ; r(t) = (x(t) , y(t) , z(t)) , 



(1) 



where the subscript t denotes the time-derivative, (x, y, z) 
are the variables of the 3D phase space, v is a given 
function of these variables, and a is a control parameter. 



Driven damped nonlinear oscillators, Duffing's equations 
for nonlinear mechanical vibrations, the Lorenz equations 
describing 2D convection in a heated fluid, nonlinear elec- 
tronic circuit equations, etc., are some well-known exam- 
ples [4] of Eq. (1). 

A phase trajectory of the dynamical system of Eq. (1) 
can be viewed as a space curve generated by the three- 
dimensional vector r(t) parametrized by the time t. From 
Eq. (1), we have |r t | = |v| = v = s t , giving s(t) = J v dt 
as the arc length on the space curve. Thus the unit tan- 
gent vector is given by T = r t /v. In accordance with 
the usual formalism [5] for a space curve, we define the 
orthogonal right-handed triad of unit vectors (T, N,B), 
where N and B denote, respectively, the normal and bi- 
nomial unit vectors on the curve. The Frenet-Serret (FS) 
equations can be written in terms of the variable t as 



T t = vnN ,N t = -MT + iiTB,B( = -DTN. 



(2) 



where the curvature k and the torsion r are functions of 
s and determine the local geometry of the trajectory. It 
can be shown that [5] 



k(s) = \r t x r tf |/ |r t | 3 , 

t(s) = r t • (r tt x r m )/\r t x r tt \ 2 . 



(3) 
(4) 



Intuitively, the curvature measures the departure of a 
curve from a straight line, while the torsion measures its 
non-planarity. 

To understand how the anholonomy of a trajectory 
arises [6], we rewrite Eqs. (2) as F t = £ x F, where 
F stands for T, N, or B, and £ = -v kB + vtT. This 
shows that as one moves on the trajectory r(t), the FS 
triad (T, N, B) rotates with angular velocities vk and 
vt around B and T respectively. In the (N, B) plane, 
we may introduce two orthogonal unit vectors u and w 
such that the triad (T, u, w) does not rotate around T. 
This is achieved by using a Fermi- Walker parallel trans- 
port of any vector P moved along the curve according to 
SP/6t = Pt — v k(B x P). Thus, as one moves from t = 
to t = T, a geometric phase <I>t = J vt dt — J^/ a \ t ds 
develops between its natural frame (N, B) and the "non- 
rotating frame" (u, w). The anholonomy $t can also 
be interpreted in another way: Representing the rota- 
tion of the triad using Euler angles (6, <p, tp) yields [6] 
$t = 27T — J sin dO J dip. For a periodic trajectory, the 
second term is just the solid angle subtended by the area 
enclosed by the closed path -f(s) traced out by the tan- 
gent indicatrix [5] on the unit sphere S 2 . The same result 
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holds good for a non-periodic trajectory as well, since it 
can always be closed using a geodesic [7] on the sphere. 

As is clear, the anholonomy $t of a trajectory also 
characterizes its non-planarity. Since both periodic as 
well as chaotic trajectories occur in general, a more useful 
quantity is the geometric phase per unit time, defined 
over a long time T, 



r>j< = 



rdt 



(5) 



where T is much larger than all the natural time scales 
in the problem. (f>T as defined above can also be thought 
of as the mean angular velocity of the FS triad. As we 
show below, it turns out to be a good order parameter 
for both periodic and chaotic trajectories. 

We now apply this formulation to a driven, damped, 
bilinear oscillator [9] characterized by two different fre- 
quencies wi,2 for positive and negative displacements re- 
spectively: 



x t t + 20x t — -oj\ 2 x + f cos ut , (0 > 0) . 



(6) 



This can be written in the form of Eq. (1) by choos- 
ing r = (x,y,z) = (x,x t ,x tt ). For = 0, / = , the 
oscillator has a frequency Wbi = 2ojiw 2 /(^i + uj 2 )- The 
system is piecewise linear, and analytic solutions can be 
obtained for x > and x < 0. The discontinuity at the 
origin makes it essentially nonlinear, a = W2/0J1 being 
the nonlinearity parameter. Without loss of generality, 
we choose the units of t and x such that w = 1 and / = 1 . 
The two key control parameters are then u)\ and uo 2l or, 
equivalently, a and oj\. The results to be presented are 
for fixed w\ — 1.5, although we have investigated the 
full two-dimensional parameter space [8]. For a linear 
damped driven oscillator {lo\ = 102), it is easy to show 
that r calculated using Eq. (4) is identically zero and 
hence the attractor is a planar limit cycle. 

Our numerical routine uses the analytical solution for 
x(t) in each half of the phase space. Starting with an ini- 
tial state x > , xt — (so that the oscillator frequency 
is lo\), we use very small time increments to determine 
precisely the time and the velocity when x reaches 0. 
These are used as the new initial conditions for dynam- 
ics in the region x < with oscillator frequency u>2, and 
the process is repeated. The availability of an analytic 
solution is particularly useful for computing the torsion 
of the phase trajectory using Eq. (4). 

We have carried out a detailed exploration [8] of the 
system in the parameter space, and calculated the local 
phase space variable x, the local geometric variable r, 
the long-time average anholonomy <j>T , as well as 4>i , the 
average anholonomy over a single period of the driving 
force. 
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FIG. 1. Bifurcations and geometrical resonances: The at- 
tractor at Poincare points x„ and the geometric phase factors 
4>i and 4>t as the nonlinearity parameter a is varied. Plots 
for a > 1 and a < 1 are shown separately to highlight the 
differences in the magnitude of the anholonomy. 

The bilinear oscillator exhibits period-doubling tran- 
sitions, followed by chaotic dynamics which is reflected 
by 4>i in addition to phase space variables. A close look 
at Fig. 1 shows that all changes in the attractor in 
phase space are reflected in (f>T as well, either as a dis- 
continuity in </>t or as a change in its slope dipr/da. For 
example, at a = 2/3 (o>2 = 1) and a w 4.25 there is a 
sudden change in the Poincare plot x n (due to coexist- 
ing attractors changing their basin of attraction) which 
is accompanied by a jump in <pT- A sequence of period- 
doublings at a w 3.4, 4.25, 4.7... as well as a resonance at 
a w 2.0, where u>m = to, are accompanied by changes in 
depr/da. 

A striking feature is the emergence of a series of 
new geometric transitions of the attractor that are not 
captured by the standard bifurcation diagram for x n . 
While x n varies smoothly, <pi and (j)T show a jump at 
a Rj 1.25, 3.0, 3.4 , reminiscent of a first-order phase tran- 
sition. In contrast, at a rs 0.61 , x n and <f>T are both 
continuous, but d(f>T/da has a jump, reminiscent of a 
second-order phase transition. These transitions are ro- 
bust, as the threshold parameter values for both first and 
second order transitions lie on a smooth one-dimensional 
curve in two-dimensional parameter space [8] . Based on 
our detailed numerical studies, there is ample justifica- 
tion for regarding <f>T as an order parameter for such ge- 
ometric transitions. 

The first order transitions result in a discontinuous 
change in ^ by u ( which is set to unity here) in only 
one of the branches of x n . As a result, we see discontinu- 
ities in (j>T by exactly 1/n, where n is the period of the 
attractor in units of 2-k/lo. Therefore, the discontinuities 
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in 4>t correspond to a full rotation of the FS triad by 27T 
in a full period of the driving force. These transitions 
can thus be viewed as a resonance effect between the av- 
erage rotation frequency of the FS triad and the driving 
frequency. Higher order resonances where the FS triad 
rotates by 27rm (m = any integer) were also frequently 
observed [8] for a < 1 . In the neighborhood of a second- 
order transition, the mean anholonomy is very large as 
the FS triad executes many rotations (not necessarily an 
integer number) during one period of the driving. Such 
sudden responses in the geometry of the trajectories that 
are not captured by usual bifurcation diagrams will be 
referred to as geometrical resonances (GR). 
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FIG. 2. Geometrical localization near first-order (left) and 
second-order(right) geometrical transitions. Time series for 
r (thick lines) and x (dashed lines; full z-scale not shown), 
shown over one driving cycle. The horizontal line corre- 
sponds to x = 0. Unlike x and x t , r is sensitive to the 
discontinuity in the acceleration at x = 0. /,F rom top to 
bottom: on the left, a — 1.1, 1.2, 1.25, 1.3; on the right, 
a = 0.66, 0.62, 0.61, 0.59. 

Further insight into these intriguing phenomena can 
be gained by studying the local dynamics of both the 
phase space and geometrical variables as shown in Fig. 
2 . As we discuss below, the key to the various effects 
is the localization characteristics underlying the torsion 
r, which we will refer to as the T-mode. We first discuss 
the dynamics near the first-order transition, given on the 
left in Fig. 2. An interesting aspect is the temporal local- 
ization of the r-mode corresponding to large geometrical 
changes in very short time intervals, as the nonlinearity 
increases from a = 1.1. At a critical parameter value, 
it becomes a ^-function and then flips, i.e., the localized 
r changes its sign, after which the r-dynamics that was 
"bi-directional" (i.e., r had both signs) becomes "uni- 
directional", with the FS triad rotating in one direction 



alone, as can be inferred from the figure. The geometrical 
transition appears to be linked to subtle changes in the 
dynamics. Here the local curvature vanishes as r t and 
y u become collinear. In terms of phase space variables, 
the force and its two derivatives are very small and the 
particle is subjected to a small constant force for a long 
time interval. 

One of the interesting aspects of a first-order GR tran- 
sition is the fact that it is often associated with a period- 
doubling bifurcation; in that case it is a "twin" transi- 
tion, where a jump in <fii before period-doubling is al- 
ways accompanied by another jump (of identical magni- 
tude but opposite sign) in </>i after the period-doubling. 
The twin aspect of the GR was seen throughout the two- 
dimensional parameter space, and may be universal. 

In contrast to the first-order case, the dynamics near a 
second-order transition corresponds to an enhancement 
in the localization length of the r-mode. as described on 
the right in Fig. 2. In addition to this derealization, the 
mode also becomes predominantly unidirectional (posi- 
tive r) due to the disappearance of the strongly localized 
(negative) r. After this stage, the mode becomes bi- 
directional again. In view of the absence of "flipping", 
4>t remains continuous, while its derivative exhibits a 
discontinuity. 
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FIG. 3. Time series for r (thick line) and x (thin line) in 
the chaotic state for a = 5. The middle and the bottom plots 
show the blowups. 

Localization of r implies that the non-planarity of the 
attractor is significant only in short intervals of time. 
This nonlinearity-induccd localization of the geometri- 
cal variables when the phase space variables continue to 
exhibit oscillatory dynamics is an interesting new phe- 
nomenon that we will refer to as geometrical localiza- 
tion. In addition to their existence near GRs, localized 
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r-modes are seen frequently as sudden, seemingly ran- 
dom bursts in the chaotic time series as clearly seen in 
Fig. 3. They are tied to very subtle changes in the dy- 
namics associated with the derivatives of x. 

We had stated earlier that 4>t is an appropriate order 
parameter even for dynamics on a chaotic attractor. This 
is borne out by the fact that <pT converges to a well- 
defined asymptotic value for very large sampling times t, 
as shown in Fig. 4. 
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FIG. 4. Convergence of <f>T for the attractor ( i.e. large 
t) for various values of a in the periodic (left) and chaotic 
(right) dynamics. The two curves in each of the plots on the 
left show 4>t before and after GRs, respectively. 



systems using the FS equations provides a new method 
which adds geometrical features to the usual phase space 
description of the complex dynamics of classically chaotic 
systems. Non-planarity of the attractor is related to an- 
holonomy, which, in turn, is shown to be an appropriate 
robust order parameter that characterizes the geometry 
of the attractors. Adding to the richness of the bifurca- 
tion diagram is the prediction of first and second-order 
transitions in the geometry of the attractor. These tran- 
sitions are related to the localization characteristics of 
the torsion. This suggests the novel concept of geomet- 
rical localization near a transition, which emerges due to 
the fact that significant deviation from planar dynamics 
exists only over small time intervals. We thus relate non- 
linearity not only to anholonomy but also to localization. 

We believe that several of the features we have seen in 
our geometric description of the bilinear oscillator may 
be generic and found in other nonlinear systems as well. 
These were seen throughout the two-dimensional param- 
eter space [8] , and also for other choices of phase space 
variables such as r = (x, Xt, t) and (x, xt, cos u)t). Our 
formulation can also be extended to D > 4. For instance, 
when D = 4, the anholonomy can be found from the ex- 
pression for the angular velocity of a tetrad that can be 
appropriately defined for the system. We hope that our 
study will stimulate a new line of research relating geom- 
etry and dynamics in nonintegrable systems in general. 

The research of IIS is supported by National Science 
Foundation Grant No. DMR 0072813. 



The average anholonomy fluctuates before converging 
to a fixed value, reminiscent of the critical fluctuations of 
an order parameter near a thermodynamic phase transi- 
tion. The fact that 4>t converges to a unique asymptotic 
limit even for chaotic trajectories establishes the resilient 
aspect of the geometric phase in the classical bilinear os- 
cillator. This tolerance to fluctuations may be a universal 
characteristic of geometric phases in both quantum and 
classical systems which has led to the exciting possibility 
of geometric quantum computation [10]. 

It is interesting to note that, by defining a complex 
vector M = (N + iB) / y/2, a short calculation using the 
FS equations (2) shows that the classical geometric phase 
3>t = J r / a \ T~ds = i J r , . M* • M. s ds. This expression has 
exactly the same form as the quantum geometric phase 
found by Berry [1], when the classical vector M is re- 
placed by a quantum state |n(R)). In view of the gen- 
eral robustness of the geometric phase, we believe that 
our formalism that uses anholonomy may provide an al- 
ternative formulation for studying the classical-quantum 
correspondence. 

In summary, our space curve formulation of dynamical 
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